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Antiferro-quadrupole state of orbital-degenerate Kondo lattice model with f 2 configuration 
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To clarify a key role of / orbitals in the emergence of antiferro-quadrupole structure in PrPb3, we inves- 
tigate the ground- state property of an orbital-degenerate Kondo lattice model by numerical diagonalization 
techniques. In PrPb3, Pr 3+ has a 4/ 2 configuration and the crystalline-electric-field ground state is a non- 
Kramers doublet T3. In a j-j coupling scheme, the T3 state is described by two local singlets, each of 
which consists of two / electrons with one in Y7 and another in Fg orbitals. Since in a cubic structure, 
T7 has localized nature, while Tg orbitals are rather itinerant, we propose the orbital-degenerate Kondo 
lattice model for an effective Hamiltonian of PrPb3. We show that an antiferro-orbital state is favored by 
the so-called double-exchange mechanism which is characteristic of multi-orbital systems. 
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It is currently one of the central issues in the research field 
of condensed-matter physics to unveil novel magnetic phases 
of strongly correlated electron systems with active orbital de- 
grees of freedom. It has been a common understanding that 
competition and interplay among spin, charge, and orbital de- 
grees of freedom cause diverse ordering phenomena involv- 
ing multiple degrees of freedom, as frequently observed in d- 
and /-electron systems. 1 ' 2 ) In the case of /-electron systems, 
spin and orbital are tightly coupled with each other due to the 
strong intra-atomic spin-orbit interaction. To describe such a 
complex spin-orbital state, the /-electron state is usually clas- 
sified in terms of multipole degrees of freedom. 

A rare-earth compound PrPb3 with a simple AuCu3-type 
cubic structure has attracted great interest as a typical material 
that exhibits antiferro-quadrupolar (AFQ) ordering. In fact, 
this compound undergoes a second-order transition at 0.4 K, 3) 
which has been confirmed to be a non-magnetic but an AFQ 
transition. 4, 5) In PrPb3, Pr 3+ has a 4/ 2 configuration, and the 
crystalline-electric-field (CEF) ground state is a non-magnetic 
non-Kramers doublet T3 with a magnetic triplet lying 19 K 
above the ground state. 6 ' 7) The T3 state carries O® and 0\ 
quadrupole moments. Thus, the low-temperature property is 
governed by quadrupole degrees of freedom. 

Regarding the H-T phase diagram of PrPb3, the so-called 
reentrant phase diagram has been obtained, in which the tran- 
sition temperature goes up with increasing the field but turns 
to decrease and the ordered phase closes at a low field. 7 ) This 
reentrant behavior has been well reproduced phenomeno- 
logically based on a mean-field theory assuming a simple 
two-sublattice ordered structure. 7) However, recent neutron 
diffraction measurements have revealed that the quadrupole 
ordered structure is modulated in space, instead of a sim- 
ple two-sublattice structure. 8 ) In principle, such a long-period 
ordered structure could emerge because of significant long- 
range quadrupole interactions, although the origin of the long- 
range interactions is not clear. 

So far, there have been no theoretical efforts to understand 
AFQ structure of PrPb3 from a microscopic viewpoint. In this 
paper, we propose an orbital-degenerate Kondo lattice model, 
which is obtained on the basis of a j-j coupling scheme, as 
an effective model for PrPb3. We investigate the ground-state 




Fig. 1. Schematic view of T% and T7 orbitals. 

property of the model by using exact-diagonalization tech- 
niques. It is found that an antiferro-orbital state emerges due 
to the so-called double-exchange mechanism which is in gen- 
eral relevant to multi-orbital systems. 

First we explain the construction of an effective model for 
PrPb3. In the j-j coupling scheme, we first include the strong 
spin-orbit interaction, and we accommodate / electrons in the 
lower sextet with the total angular momentum j=5/2. Under 
the cubic CEF effect, the sextet is split into a T7 doublet and 
a Ts quartet. To distinguish two Kramers doublets in the Tg 
quartet, it is useful to introduce two orbitals, while spin is also 
introduced to represent two states in each Kramers doublet. 
Note that the Tj doublet gives another orbital. The schematic 
views of Tg and orbitals are depicted in Fig. 1. Since we 
accommodate / electrons in the level scheme of the one /- 
electron state, we refer to the level scheme of CePb3, which 
is a 4/ 1 compound with the same lattice structure with that of 
PrPb3. In CePb3, it has been found that is the ground state 
and Ts is the excited state. 9 ) Thus, for PrPb3, we accommo- 
date two / electrons in this level scheme. 

In PrPb3, the CEF ground state is the non-Kramers doublet 
T3. In the j-j coupling scheme, the Ts state is described by 
two local singlets, each of which is composed of two electrons 
with one in Tj and another in Tg orbitals. Here we note that 
Ts orbitals carry quadrupole degrees of freedom. Taking ac- 
count of the formation of local singlets, an antiferromagnetic 
(AFM) coupling should be effective between electrons in 
and Ts orbitals, although the Hund's rule coupling causes a 
ferromagnetic (FM) coupling. Thus, here we introduce the 
AFM coupling as an effective interaction to involve a high- 
order CEF effect B§ which can not be included in the .7=5/2 
Hilbert space in the j-j coupling scheme. 
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Concerning the itinerancy and localized nature of orbitals, 
we consider an /-electron hopping through the sigma bond. 
In a cubic structure, due to the spatially anisotropic shape of 
orbital, orbital is localized, while Tg orbital is itinerant in 
the xy plane and Tg orbital is itinerant in all three directions. 
Thus, to consider an effective model, we assume that electron 
in T-i orbital is localized, leading to a localized spin. Note that 
T 8 electron is itinerant and couples with localized T 7 spin due 
to the effective AFM interaction, which can be regarded as an 
analog of the Kondo coupling. 

Taking into account these situations, we obtain an orbital- 
degenerate Kondo lattice model as an effective Hamiltonian 
for PrPb3, given by 

t,t'ct i 

+ U ^2 P'^Pirl + U'^Y] PiaPi(3 
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where f\ ra is the annihilation operator for electron with 
spin cr(=t, j) in orbital r(=a,(3) at site i, p ira =fl a f ir(T , 

Pir=Y,a Pirai Sir 8 = (l/2) Ysaa'T /L^'/^' where ^ 

are Pauli matrices, and Sir 7 is the spin- 1/2 operator for I\ 
spin. The summation of (ij) is taken for nearest neigh- 
bor sites in the cubic lattice. The hopping amplitudes are 
evaluated from the overlap integral between /-orbital wave- 
functions in adjacent sites, which are given by t* a =3£/4, 
t^=t^ a =-V3t/4, t%p=t/4 for the x direction, tl a =3t/4, 
t^ /3 =t^ a =v / 3t/4, t*p=t/4 for the y direction, and tpp=t, 
t z aa =t z a(3 =t z p a =0 for the z direction, where t=(3/7)(f fa). 
Hereafter, t is taken as the energy unit. In the second term, 
Jk is the Kondo coupling between I\ spin and T$ electron. 
The rest terms are interactions among T$ electrons: U, U' , 
J, and J' denote intra-orbital, inter-orbital, exchange, and 
pair-hopping interactions, respectively. Note that the relation 
U =U'+J+J' holds, which originates from the rotational in- 
variance in the orbital space, and J=J' is assumed. 2 ^ 

We analyze the model (1) by numerical diagonalization. 
Since the size of the Hilbert space becomes so large as 32^ 
due to orbital degree of freedom, where N is the number of 
sites, it is rather difficult to enlarge the system size. However, 
the method is advantageous to grasp the ground- state prop- 
erty, such as orbital structure, in an unbiased manner. In the 
present work, first we study a 2x2 square four-site system 
in the xy plane. Then, taking account of the characteristics 
grasped within the four-site system, we proceed to a 2 x 2 x 2 
cubic eight- site system. In this paper, we set U' /W=2, where 
W is the band width, and investigate the dependence on Jk 
and J. Note that the band width is W=A for the square system 
and W=6 for the cubic system. 

First, we show the results for the four-site system. The main 
result is summarized in Fig. 2(a), which is the ground-state 
phase diagram in the ( Jk, J) plane. There three types of com- 
peting spin-orbital configurations are observed. When Jk and 
J are small, we find a ferro-orbital (FO) state with an AFM 
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Fig. 2. Four-site results, (a) Ground-state phase diagram in the (Jk,J) 
plane. Inset denotes schematic view of electron configuration in each 
phase, (b) T(q) as a function of J at Jk=0.6. (c) FO and AFO struc- 
tures, (d) S(q) as a function of J at Jk=0.6. (e) C s as a function of Jk 
at J=0. Note that C s (i) takes the equivalent value at every site due to the 
translational symmetry. 

configuration in each of T 7 and T 8 orbitals, while spins in 
and T 8 orbitals are antiparallel at each site due to Jk. With 
increasing J, an antiferro-orbital (AFO) state occurs, and the 
spin state turns to be FM at a larger J than the transition point 
from FO to AFO. It should be noted that even though the spin 
state is characterized by AFM or FM in each of T7 and Tg 
orbitals, a local singlet is formed due to Jk and the ground 
state is totally non-magnetic. 

Let us here discuss the orbital state. In order to determine 
the orbital structure, it is useful to introduce new operators 
for £ and r\ orbitals, which are given by linear combinations 
of the original operators, such as 

fi iG =cos(0i/2)fip a + sin(0i/2)./w, 

= - sin(0i/2)/ i/3a + cos(^/2)/ iaa , (2) 

where 0\ characterizes the orbital shape at each site. The opti- 
mal {#i} is determined so as to maximize the orbital structure 
factor, defined by 

T(q) = ^(^T k V q - (j - k) /^, (3) 
where f.^MU^-fU^)/ 2 and (•••) denotes the 
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expectation value. In Fig. 2(b), we show T(q) as a function of 
J at Jk=0.6 with 0\=9. For small J, the dominant component 
is T(0, 0) with 0=tt, indicating FO state. In the FO state of the 
present square system in the xy plane, Tg orbitals are favor- 
ably occupied to gain kinetic energy, since orbital extends 
to adjacent sites, as shown in Fig. 2(c). With increasing J, the 
dominant component changes to T(7r,7r) with 0=tt/2, indi- 
cating AFO state in Fig. 2(c). We note that the present square 
lattice structure is reflected in the orbital shape. In this con- 
text, it is crucial to proceed to the eight-site system to discuss 
the orbital structure in the cubic lattice. 

Concerning the spin state, we measure the spin structure 
factor of T 8 electrons, defined by 

^(q) = E^s^r 8 )e iq - (j - k) /^. (4) 

As shown in Fig. 2(d), with increasing J, 5(0, 0) is increased, 
while 5(7r, tt) is reduced, so that the dominant spin correla- 
tion changes from AFM to FM. Here we note that the motion 
of T 8 electrons leads to a FM spin arrangement. Namely, the 
double-exchange mechanism is effective, which is character- 
istics of multi-orbital systems. 

We also measure the on-site spin correlation, defined by 

C 8 (i) = (S? rr S? r8 ). (5) 

In Fig. 2(e), the Jk dependence of C s at J=0 is shown. It 
is obvious that C s =0 at Jk=0, since there is no correlation 
between I\ spin and electron. With increasing Jk, C s de- 
creases and gradually approaches —1/4, indicating the stabi- 
lization of the local singlet. It is found that even when J is 
increased, C s keeps a value near —1/4 (not shown), indicat- 
ing the robust formation of the local singlet. Here it is worth 
noting that as Jk increases and the local singlet is stabilized, 
the AFO phase tends to extend to the region of small J, as 
shown in Fig. 2(a). Namely, the double-exchange mechanism 
becomes significant due to Jk. 

Now we move on to the results for the eight-site system. In 
Fig. 3(a), we show the 6 dependence of T(q) at Jk=1 and 
J=0. We find that T(tt, tt, tt) is dominant, indicating AFO 
state. Concerning the orbital shape, it is observed that the 
magnitude of T(tt, tt, tt) does not depend on 6. We can not 
determine the actual orbital shape, but the AFO structure with 
any 6 is possible to realize. The AFO structure with 0=0 is 
shown in Fig. 3(b). As for the spin state, S(tt, tt, tt) is found 
to be dominant. We note that the present parameter set ( Jk=1 
and J=0) is corresponding to the FO phase in the four-site 
system, as shown in Fig. 2(a). It is naively expected that the 
AFO phase extends to a broad area in the phase diagram even 
when we consider the cubic lattice. However, at Jk=1 and 
J=4, T(q) has dominant components T(tt, tt, 0) with 0=tt/2, 
T(tt, 0, tt) with #=57r/6, andT(0, tt, tt) with 0=tt/6, as shown 
in Fig. 3(c), while 5(0, 0, 0) is dominant. The (tt, tt, 0) orbital 
structure is depicted in Fig. 3(d). When J is further increased, 
the orbital structure is considered to turn to be AFO due to 
the double-exchange mechanism, suggesting a rich phase di- 
agram including competing orbital states. 

Finally, we briefly discuss possible relevance of the present 
results to the AFQ structure of PrPb3. We have shown that 
for a one-dimensional j-j coupling model with an J 1 con- 
figuration, 10) an incommensurate orbital state appears due to 
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Fig. 3. Eight-site results, (a) T(q) as a function of and (b) orbital struc- 
ture at Jk=1 and J=0. (c) T(q) as a function of 6 and (d) orbital structure 
at Jk=1 and J=4. 

the competition between itinerant and localized orbitals. By 
analogy, we expect that the competition among plural orbital 
states with different nature could cause a modulated orbital 
structure as observed in PrPb3. 

In summary, we have investigated the ground- state prop- 
erty of the orbital-degenerate Kondo lattice model to under- 
stand the quadrupole structure in PrPb3 from a microscopic 
viewpoint. We have observed several types of competing spin- 
orbital states. In particular, it has been emphasized that the 
AFO state emerges due to the double-exchange mechanism. 
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